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^ Abstract 



odd number of dimensions. We find no evidence of any instability in five or seven dimensions, but 
f~| in nine dimensions, for sufficiently rapid rotation, we find perturbations that grow exponentially in 

time. The onset of instability is associated with the appearance of time-independent perturbations 
Oh which generically break all but one of the rotational symmetries. This is interpreted as evidence 

^ for the existence of a new 70-parameter family of black hole solutions with only a single rotational 

I symmetry. We also present results for the Gregory-Laflamme instability of rotating black strings, 

demonstrating that rotation makes black strings more unstable. 
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1 Introduction 

In four spacetime dimensions, the Kerr solution is believed to be stable against small perturbations of 
the metric: if it is perturbed then it will eventually settle down to another Kerr solution with slightly 
different mass and angular momentum. Although this has not been proved rigorously, evidence in 
favour of this picture arises from studies of linearized gravitational perturbations of the Kerr solution 
|35[ [37] , and from numerical studies of black hole formation [3l| . 

In higher dimensions, the stability of the Schwarzschild solution against linearized gravitational 
perturbations has been established [I H 122 ^ [20]. However, it has been argued that the D-dimensional 
generalization of the Kerr solution, discovered by Myers and Perry [30j, should exhibit an instability 
for sufficiently large angular momentum. The MP solution is parameterized by its mass M and 
\_{D — l)/2j angular momenta (where [J stands for the smallest integer part). In Ref. [8], it was 
shown that if only one angular momentum is non-vanishing then, in the limit in which this angular 
momentum is very large (using a scale set by M), the black hole locally resembles a black brane. 
However, black branes suffer from the classical Gregory-Laflamme instability [12]. This suggests that 
"ultraspinning" MP black holes will also be unstable. The argument can be generalized to cases in 
which several angular momenta become large while others remain small [8]. 

To confirm the existence of the ultraspinning instability, and to identify the critical angular mo- 
mentum beyond which an instability exists, requires a linear stability analysis of the Myers-Perry 
solution. For the Kerr solution, the equations governing linearized gravitational perturbations can be 
decoupled and reduced to a single scalar wave equation, which renders a stability analysis tractable 
|35] . However, no analogous decoupling has been discovered for the MP solution so the equations to 
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be studied consist of many coupled linear partial differential equations, making the problem seem very 
difficult. 

Recently there was significant progress with this problem. Consider a MP solution with a single 
non-vanishing angular momentum J. For given M, the ultraspinning instability is expected to occur 
only when J exceeds some critical value Jcrit determined by M. A solution with J = Jcrit is expected 
to admit a stationary zero-mode, i.e. a time-independent perturbation that is regular on the future 
horizon and vanishes at infinity. Ref. [1] determined Jcrit and the stationary zero-mode for D = 7,8, 9. 
This was done by writing down a certain Ansatz for the form of the metric perturbation and solving 
the resulting PDEs numerically using a novel (in this context) method that we shall explain below. 

We regard the work of Ref. [3] as very strong evidence for the existence of an ultraspinning 
instability beyond J = Jcrit- However, no actual instability, i.e. a perturbation growing in time, was 
demonstrated. In this paper, we shall demonstrate that some MP black holes do admit perturbations 
that grow exponentially in time, thereby providing the first example of a linearized gravitational 
instability of an asymptotically flat vacuum black hole solution. 

We shall exploit the idea introduced in Ref. [24J of considering MP solutions with enhanced symme- 
try. The generic MP solution has isometry group Rx [7(1 )" where R corresponds to time translations 
and n = [{D — l)/2j. However, this is enhanced when some of the angular momenta coincide. In 
particular, for odd D, the MP solution with all angular momenta equal has a much larger R x U{N-\-l) 
isometry group, where D = 2N + 3. Furthermore, the solution is cohomogeneity-1, i.e. it depends 
only on a single coordinate. The metric involves a fibration over complex projective space CP^ . 
Gravitational perturbations of this solution can be decomposed into scalar, vector and tensor types 
according to how they transform under isometrics of CP^ . The tensors, which exist only for N > 2 
{D > 7), were studied in Ref. [24j and no evidence of any instability was found. The special case of 
D = 5, for which only scalar perturbations exist, was studied in Ref. [29]. Again, no evidence of any 
instability was found. 

In this paper, we shall study scalar-type perturbations of these cohomogeneity-1 MP black holes. 
The symmetries enable the problem to be reduced to coupled linear ordinary differential equations 
(ODEs) which we solve numerically. We find no evidence of any instability for D = 5 (consistent with 
Ref. [29]) or D = 7. However, for D = 9, when J exceeds a certain critical value Jcrit, there is a 
perturbation that grows exponentially in time, i.e. an instability. We believe that such an instability 
will exist for all (odd) D > 9 although we have demonstrated this only for D = 9. 

As expected, the onset of instability is associated with the appearance of a stationary zero-mode 
of the MP solution with J = Jcrit- This zero-mode is interesting for another reason. It has been 
proved that a stationary, rotating black hole must admit a rotational isometry (i.e. a U{1) isometry) 
|17| I26j. However, all known higher-dimensional black holes have multiple rotational isometrics (e.g. 
MP black holes have [{D — l)/2j commuting U{1) isometrics), i.e. more symmetry than one expects 
on the basis of general arguments. Therefore, it has been conjectured that there exist solutions with 
less symmetry than any known solution, specifically solutions with a single rotational symmetry [32]. 

It was proposed in Ref. [32j that one could seek evidence for the existence of such solutions in the 
same way that the first evidence was obtained for the existence of non-uniform black string solutions. 
For black strings, the static zero-mode associated with the onset of the Gregory-Laflamme instability 
was conjectured to describe the "branching off' of a new family of non-uniform black string solutions 
from the already known branch of uniform solutions [18] . Perturbative [H] and numerical p8j work 
subsequently confirmed that this was correct. For rotating black holes, the idea proposed in Ref. |32] 
is to look for a stationary zero-mode of a MP solution. By analogy with the black string example, 
this could be interpreted as the branching off of a new family of solutions. If the zero-mode preserves 
only a single rotational symmetry then this would be evidence for the existence of new black holes 
with just one rotational symmetry. 

In our case, the stationary zero-mode generically preserves only a single rotational symmetry. 
Therefore we conjecture that there exists a family of stationary black hole solutions with just a single 
rotational symmetry, that bifurcates from the cohomogeneity-1 MP black hole solution at J = Jcrit- 
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In fact, we find not just one stationary zero-mode, but a large family, corresponding to all scalar 
harmonics on CP^ with a certain eigenvalue. For = 9, we shall argue that the new family of 
solutions will involve 70 independent parameters, considerably more than the 5 parameters required 
to specify the MP solution! If correct, this implies that any hope of specifying higher-dimensional 
black holes uniquely using just a few parameters is bound to fail. Of course, the new black hole 
solutions may turn out to be unstable themselves. 

Recently, Ref. [7j has provided other evidence for the existence of higher-dimensional black holes 
with a single rotational symmetry. Using the "blackfold" approach of Refs [Sj |6] , approximate solutions 
were constructed for D > 5 that describe "helical" black rings. The blackfold approximation is based 
on the fact that higher-dimensional black holes can have widely separated horizon scales. The results 
of the present paper concern black holes that lie outside the regime of validity of this approximation. 
Our results complement those of Ref. [7] because we are presenting evidence for topologically spherical 
black holes with a single rotational symmetry whereas Ref. [7] considered black rings. 

A difference between our results and the expectations of Ref. [8] concerns the nature of the sta- 
tionary zero-mode. In our case, this perturbation breaks all but one rotational symmetry. However, 
in Ref. [8], it is argued that for single-spinning black holes the stationary zero- mode should preserve 
the symmetries of the background. The results of Ref. |1] confirmed this expectation. Singly spinning 
black holes exhibit symmetry enhancement for D > 6. They are cohomogeneity-2 with isometry group 
R X U{1) X SO{D - 3) where SO{D - 3) has S^'^ orbitsQ The zero-mode constructed in Ref. [4j 
preserves all of the symmetries. Our results raise the question of whether singly spinning MP black 
holes might admit further stationary zero-modes that break some of their symmetry, and provide 
further evidence for new black hole solutions with reduced symmetry. 

A corollary of our approach is the first data for the Gregory-Laflamme instability [12i| for rotating 
black strings. We consider black strings obtained as the product of a cohomogeneity-1 MP black hole 
with a flat direction. We argue that such solutions are always classically unstable. Our numerical 
results demonstrate that the string becomes more unstable as the angular momentum increases: the 
instability becomes stronger (i.e. it occurs on a shorter time scale) and the critical wavelength of 
unstable modes decreases, as the angular momentum increases. 

An important feature of cohomogeneity-1 MP black holes is that they exhibit an upper bound 
on their angular momentum. Solutions saturating this bound are extreme black holes. Because of 
this upper bound, it is not obvious that such black holes should exhibit the ultraspinning behaviour 
discussed in Ref. 0. Therefore we should explain why we undertook this project. 

Ref. [8] observed that the transition to "black brane-like" behaviour of singly spinning MP black 
holes can be seen in the thermodynamics: a certain thermodynamic quantity changes from "Kerr black 
hole-like" to "black brane-like" at a critical value of the spin. This was generalized in Ref. [Ij , where 
it is argued that the quantity of interest is the Hessian matrix Hij = {d'^{—S)/dJidJj)M where S is 
the black hole entropy. This is positive deflnite for small angular momenta. However, it can become 
indefinite as the angular momenta increase. Hence, one can define an ultraspinning black hole to be 
one for which this Hessian fails to be positive definite. 

For D = 5 MP black holes, Hij is always positive definite [3] so such black holes are never ultra- 
spinning. However, for D > 6, singly spinning MP black holes with large enough angular momentum 
are, of course, ultraspinning and the numerical results of Ref. [3] supply strong evidence that the 
instability appears only when Hij fails to be positive definite. 

The present paper was motivated by the observation that, for D > 7, cohomogeneity-1 MP black 
holes do satisfy the ultraspinning criterion once J exceeds a critical value Juitra (for given M). Hence, 
there might be an ultraspinning instability for J > Jcrit where Juitra < Jcrit < ^extreme- Our results 
show that there is no instability for D = 7 but an instability does occur for D = 9 and, we believe, 

^ One can decompose metric perturbations of these solutions into scalar, vector and tensor types using this SO{D — 3) 
symmetry. The tensors, which exist only for D > 7 have been studied previously in Ref. 23 and show no evidence of any 
instability. However, tensor perturbations arise from deformations of the S^~* part of the metric, whereas the expected 
ultraspinning instability should arise from perturbations of the metric transverse to S*^"*. The results of Ref. [4] indicate 
that the instability should be a scalar-type perturbation. 
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for (odd) D>9. 

This paper is organized as follows. Section [2] describes the cohomogeneity-1 black holes that we 
shall study. In Section [3j we give a brief discussion of our approach and present our results. In 
Section [4j we show that thermodynamics can be used to explain many of our results. We also explain 
the relevance of the Hessian Hij and present plots exhibiting the region of parameter space for which 
MP black holes are ultraspinning according to the definition above. The technical details of our work 
are presented in the later Sections [5] and [6j and in the Appendices. 

Note added. As the present paper was nearing completion, Ref. |33j appeared. This paper 
performs a fully nonlinear numerical evolution of a perturbed D = 5 singly spinning MP solution. An 
instability was found for sufficiently large spin. It is not clear how this very interesting work relates to 
the instability under discussion here since D = 5 MP solutions are never ultraspinning (according to 
the definition above). Furthermore, the perturbation considered in Ref. [33] breaks the symmetry in 
the direction of rotation of the black hole whereas we (and Ref. 0]) consider only perturbations that 
preserve this symmetry. We suspect that Ref. |33| has found an instability of D = 5 MP solutions of 
a qualitatively different form from the one discussed in this paper. It would be interesting to know 
whether a similar instability occurs for D > 5. 



2 Cohomogeneity-1 Myers-Perry black holes 

The Kerr solution was extended to higher dimensions by Myers and Perry ^30j . The Myers-Perry family 
can be parameterized by a mass radius parameter and [{D — l)/2\ angular momentum parameters 
Oj. In the particular case of equal angular momenta, = a, the solution in odd dimensions D = 2A^-|-3 
is cohomogeneity-1. The metric can be written asj^ 

ds^ = -f{rfdt^ + g{rfdr^ + hirfid^j + Aadx" - n{r)dtf + r'^gabdx'^dx^ , (1) 



' \^ J.2N ^ ^2N+2 J ' ) ' J.2N+2 

2N 



g[r)h{r) ' r'^^h'^ ' 

where Qah is the Fubini-Study metric on CP^ with Ricci tensor Rab = 2(A^ + V)gab , and A = Aadx"" 
is related to the Kahler form J by dA = 2J. Surfaces of constant t and r have the geometry of a 
homogeneously squashed S"^^^^, written as an fibred over CP^ . The fibre is parameterized by the 
coordinate which has period 27r. Explicit expressions for the metric gab and Kahler potential A of 
CP^ can be obtained through the iterative Fubini-Study construction summarized in Appendix [Bj 

The spacetime metric satisfies Rfj^u = and the solution is asymptotically fiat. The event horizon 
is located at r = r+ (the largest real root of g~'^) and it is a Killing horizon of = dt + ^nd^ , where 
the angular velocity of the horizon is given by: 

^2N 

^^■f^ - 27V+2 , 2iV„2- ^"^^ 



The mass M and angular momentum J, defined with respect to d^, are [H 

M=^rlf(iV+l), J=^i^ + ^)rlfa, (3) 

where A2N+1 is the area of a unit (2A^ -|- l)-sphere. 

^The radial coordinate used here can be related to the standard Boyer-Lindquist radial coordinate of [3D] through 
-s^ + a^. 
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There is an extremality bound on the angular momentum which can be expressed as 



The solution saturating this bound has a regular, but degenerate, horizon. For fixed r+, or tm, 
ultraspinning behaviour occurs for 



Note that the range oi < a < Ogxt (for fixed rjv/) for which the black hole is ultraspinning becomes 
larger as N increases, and that ai = Cext if = 1, so there is no ultraspinning behaviour for D = 5. 



3 Strategy and Results 
3.1 Strategy 

Ref. [28] introduced new numerical techniques for determining negative modes of rotating black holes. 
In Ref. [1], these techniques were exploited to construct the stationary zero-mode expected to indicate 
the onset of an ultraspinning instability of a singly rotating MP black hole. In the present paper, we 
shall determine the stationary zero-mode indicating the onset of instability for cohomogeneity-1 black 
holes. However, our main achievement is to generalize these methods to demonstrate the existence of 
perturbations that grow exponentially in time. 

In order to explain our approach (and that of Ref. it is useful to discuss black strings. We 
shall be interested in a rotating uniform black string solution 

^^•Sstring = g^iudx^^dx" + dz'^ , (6) 

where gf^i, is a cohomogeneity-1 MP metric. Consider a transverse and traceless (TT) perturbation of 
the string of the form 

'^■^string ^ '^'^string ~^ ^ h^i,(x)dx^dx , (7) 

The TT condition is equivalent to h^i, being TT with respect to the MP metric: 

/i'^^ = , V^/i^, = . (8) 

The linearized Einstein equation reduces to 

(Al/i)^, = -VPVph^u - 2Rpp,„hP'' = -k^h^u, (9) 

where A^, is the Lichnerowicz operator for the MP background. Hence perturbations with non-zero k 
correspond to negative modes of Aj;,. The boundary conditions are that /i^^ should be regular on the 
future horizon and vanishing at infinity. 

Our strategy (following Refs. [Mill!) for studying perturbations of the black hole will be to seek 
a solution of ([9]), i.e. a negative mode of the black hole, and then vary the spin of the black hole 
until k vanishes, i.e. the negative mode becomes a zero-mode. This strategy is motivated by the 
availability of numerical techniques for solving eigenvalue equations of the form ([9|. Solutions with 
non-zero k correspond to perturbations of the black string. Therefore our method will yield results for 
the Gregory-Laflamme instability of rotating black strings as well as enabling us to search for black 
hole instabilities. 

We can Fourier analyze our perturbation in the time direction, and also decompose into Fourier 
modes in the '(/'-direction, i.e. we assume that the dependence on t and ip is given by 

V oc e-''^*+™'^, (10) 
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where m is an integer. As we shall explain in detail below, we can also decompose the perturbation 
into harmonics on CP^ . These can be of scalar, vector or tensor type. The tensors were considered 
in Ref. |24j. We shall restrict our attention to perturbations of scalar-type, which can be expanded 
in terms of scalar harmonics on CP^ . As usual, harmonics with different eigenvalue decouple from 
each other. The equations satisfied by /i^jy depend only on the eigenvalue of the harmonic in question. 
Eigenvalues of the scalar Laplacian on CP^ are labelled by a non-negative integer k (see Section [5]). 
Hence our perturbation is labelled by {uj,m,K). 

Consider the (Lorentzian) negative mode equation Q. As we have explained above, this arises 
from classical perturbations of a rotating black string. The usual approach to this problem is to 
fix (k, m, k) and to determine u. However, our approach, generalizing that of Refs. |28l 1^ will be 
to fix {u},m,K) and determine the possible eigenvalue(s) —k'^. In other words, we are determining 
the wavenumber k for which black string perturbations with given m and k have time-dependence 
associated with the given uj. We shall fix the overall scale rj\,f = 1 and determine the eigenvalue(s) 
— fc^ for fixed (w, m, k) as a increases from to extremality. If k vanishes for some value of a then 
the associated black hole admits a zero-mode with the given values of {u},m,K) (of course it must be 
checked that this is not pure gauge). 

In searching for an instability, we are looking for solutions of ^ with lm{uj) > 0. A problem 
with our approach is that we expect unstable modes to have complex lo in general, with the real and 
imaginary parts of oj related in some way. In other words, for given m, k and a, the complex quantity 
u will be a function of the real quantity k and hence the real and imaginary parts of uj cannot be 
independent. If we try to follow the above strategy for a randomly chosen complex value of uj, then 
this will not satisfy the required relation between its real and imaginary parts, and therefore our 
numerical method will not output a real value of A:. In order to locate where k vanishes we would have 
to scan over both a and, say, the real part of uj. It would be difficult to do this with high accuracy. 

We shall circumvent this problem by restricting attention to modes with m = 0, i.e. modes 
preserving the rotational symmetry of the black hole that follows from the rigidity theorem. As we 
shall explain below, there are reasons to expect that unstable modes with m = will have purely 
imaginary uj: uj = iT, F > 0, and this is confirmed by our results. Prom the black string perspective, 
for given a, if the string is unstable then we expect it to be unstable for perturbations with a range 
of k, hence there will be a range of k for which there exist solutions with F > 0. Conversely, in our 
approach, there will be a range of values of F > for which there exists a solution for k. 

In summary, we shall set m = 0, r^j = 1 and, for given (F, At, a) we shall determine the possible 
eigenvalues —k'^. Then we vary a until the eigenvalue vanishes. We have then found a black hole that 
admits an unstable zero-mode with the given values of F and k. 

3.2 Results for D = 5 

Our expectation is that, for small a, the black hole will be classically stable but the associated black 
string will suffer a Gregory-Laflamme instability [12] . Therefore, for a range of F, there should exist 
real solutions for k (corresponding to unstable perturbations of the string) but k will never vanish for 
F > 0, so the black hole is stable. For a static string, the Gregory-Lafiamme instability is an s-wave 
perturbation, which for us translates into a k = perturbation. Hence, for a rotating string, it is 
natural to expect this instability also to have k = 0. 

This is indeed what we find. The left plot in Fig. [l] shows our result for —k'^ for given a and F, 
with K, = 0. The plot with F = corresponds to a stationary perturbation of the black string. This is 
the "threshold unstable mode" at the critical wavelength beyond which the black string is unstable. 
Note that the curves exist for all values of a, i.e. the Gregory-Lafiamme instability is always present, 
it does not "switch off' as a increases. None of the curves extends to A; = so there is no indication 
of any black hole instability. 

In the right plot of Fig. [T| we give a more familiar plot of F against k for different values of a. For 
each value of a, we have a curve that takes the usual Gregory-Laflamme form. The maximum value 
of F is 10 — 20% of rjv/ and increases with increasing a. Furthermore, the range of k for which there 
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Figure 1: Results in D = 5, k = 0. We represent this mode for fixed values of Ttm (first graph) and 
0-1 I'M (second graph). This unstable mode of the string corresponds to well-known Gregory-Laflamme 
mode. 



exists an instability increases, i.e. the instability persists down to shorter wavelengths as a increases. 
Hence rotation makes the string more unstable. As usual, F — )• as fc — ?■ and as A; — )• fcc > 0. 
The mode with F = /c = has the usual interpretation of a gauge mode [13]. The mode with F = 
and A: = /cc is the threshold unstable mode associated to the onset of instability. This marks the 
bifurcation of a new family of non-uniform rotating black string solutions. These are the solutions 
constructed in Ref. [21]. They preserve the symmetries of the cohomogeneity-1 MP black hole but 
break the translational symmetry along the string. 

Note that the slope V jk approaches a common limiting value as A; — )• 0, independently of the value 
of a. This long-wavelength limiting behaviour is captured by the blackfold approach: it follows from 
Ref. [6] that F/Zc — >• l/\/D — 2 as A; — t- 0. This is consistent with our numerical results. 

We find no solution of ([9]) with k = 1. Therefore our results are consistent with stability of these 
black holes, in agreement with the results of Ref. [29]. 



3.3 Results for D = 7 

For K = 0, we have the Gregory-Laflamme instability shown in Fig. [2j This is qualitatively the same as 
for D = 5. Once again, rotation makes the string more unstable, the behaviour as A; — )■ is consistent 
with T/k ^ l/\/D — 2 independently of a, and there is a threshold mode at a critical value of k at 
which we expect a bifurcation of a new family of non-uniform strings analogous to those constructed 
in Ref. [H]. 

A new feature of = 7 is the existence of an ultraspinning regime. This occurs for o > ai, where 
ai was defined in equation ([s]). Just as in Ref. |4j, we find that a new stationary (F = 0) negative 
mode of the black hole appears at this point. This negative mode has k = 1. Our numerical results 
are shown in Fig. [3} These results demonstrate that the black string has an instability in the k = 1 
sector when a > ai. This is a new Gregory-Laflamme instability of the black string, distinct from 
the instability in the k = sector. The plots of F against k have the same qualitative shape as for 
the At = instability except that the slopes of the curves appear to vanish (for all o) as A; — )• oj^ 
Once again there is a threshold unstable mode at a critical value of k. Presumably this corresponds 
to a bifurcation to a new family of non- uniform black string solutions. In addition to breaking the 

^It would be interesting to investigate whether this behaviour can be explained using blackfold methods. 
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Figure 2: Results in D = 7, k = 0. The graphs are entirely analogous to the ones in Fig. [T| 

symmetry along the string, this mode also breaks some of the symmetry of the black hole (typically 
down to that of a generic MP black holeQ so this new family has less symmetry than the non-uniform 
strings associated to the threshold unstable mode with k = 0. 

Note that the k = 1 black string instability coexists with the k = instability. The latter is clearly 
dominant since it has much larger F and the instability exists for a larger range of k, i.e. down to 
shorter wavelengths. 

It is important to note that there is no evidence of any instability of the black hole: none of the 
curves with non-zero F extends to k = 0. In the limit A; — )• 0, solutions with non-zero F approach a 
pure gauge mode, just as for k = 0. On the other hand, the solution with F = 0, does not approach 
a gauge mode as /c — t- 0. Instead, as anticipated in Ref. [4j, it corresponds simply to a variation of 
parameters within the MP family of solutions. We shall explain below why this must be the case, and 
why we find such a mode only at a = ai. 

For D = 7, since the black hole is ultraspinning for a > ai, there is the possibility of an ultraspin- 
ning instability appearing at a = 02 > oi. However, we find no solution of equation Q for k = 2 so 
our results are consistent with stability of D = 7 cohomogeneity-1 MP black holes. 

3.4 Results for D = 9: black hole instability 

For K = 0, we have the expected Gregory-Laflamme instability of the black string. For k = 1, we find, 
as for D = 7, a new Gregory-Laflamme instability of the black string that appears at a = ai. This is 
shown in figure Fig. |4} 

The new feature that appears for = 9 is an instability in the k = 2 sector, which appears at 
a = a2 > ai. This is shown in Fig. [5] The left plot shows a new stationary (F = 0) negative mode which 
emerges from a zero-mode at a = 02- We shall prove below that this zero mode cannot correspond to 



a variation of parameters within the MP family of solutions. Furthermore, in Appendix E.2, we show 
that it cannot be a gauge mode. 

For a > 02, there is a new instability of the black string, corresponding to the curves with F > 
in the plot. But there is a qualitative difference between the left plot of Fig. [5] and our previous plots: 
the curves with F > now intersect A; = 0, i.e. we have found perturbations of the black hole that 
grow exponentially in time, that is, a classical instability of black holes with a > 02. This is our main 
result. 

* This is because k = 1 harmonics are in one to one correspondence with KiUing vector fields of CP'^ : see below. 
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Figure 3: Results in D = 7, k = 1. We represent this mode for fixed values of Ttm (first graph) 
and a/rM (second graph). It corresponds to a new Gregory-Laflamme instability of the rotating black 
string, appearing for the numerical value (oi/rjv/)"^™ = 0.5949. In the first graph, the vertical line to 
the left corresponds to the analytical prediction of ([s]), (ai/rjvf ) = 0.5946 , and the interrupted vertical 
line to the right corresponds to extremality. The second graph indicates clearly that the instability of 
the black brane does not extend to an instability of the black hole (fc = 0). 

The onset of instability is indicated by the stationary zero- mode (F = 0, = 0) at a = 02. This is 
analogous to the mode constructed for singly spinning black holes in Ref. Our main achievement 
is to demonstrate, for the first time, the existence of modes which grow exponentially with time when 
a > 02- 

The right plot of Fig. [5] shows a clear difference from our previous plots. Unlike the GL instability, 
we find that F is maximized at A; = rather than vanishing there. Hence, for the black string, the 
most unstable k = 2 modes are those with = 0, i.e. those corresponding to the black hole instability. 
For larger the black string instability "switches off' in the same way as the GL instability, with a 
threshold mode at k = kc indicating a new family of non- uniform black string solutions. 

Fig. [6] presents our result for the instability time scale of the black hole as a function of its spin. For 
o > 02, we find that F increases monotonically with a, so extreme black holes are the most unstable. 
The maximum value of F is just a few per cent of rM, so the instability is slow compared with the GL 
instability of the black string. 

Finally, we find no solutions of Q with k = 3. 

3.5 Rotational symmetries of higher-dimensional black holes 

As explained in the Introduction, the study of perturbations of higher-dimensional black holes can be 
used to investigate the possible existence of new families of black hole solutions with less symmetry 
than the known solutions. The idea proposed in Ref. |32j is to look for a stationary zero- mode of the 
black hole, which is interpreted as indicating the existence of a family of solutions branching off from 
the known solutions. 

In our case, the stationary zero- modes with k = 1, are uninteresting since (as explained below) 
these correspond to variations within the MP family. However, for D = 9, we found a stationary 
zero- mode with n = 2 that appears at a = 02, the critical value of a beyond which the black hole is 
unstable. Therefore, we have found evidence for a new family of black hole solutions that bifurcates 
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Figure 4: Results in D = 9, k = 1. The graphs are entirely analogous to the ones in Fig. [3| 



from the MP family at this point 

How much symmetry do these new solutions have? This can be inferred from the symmetry of 
the K = 2 harmonics on CP^ . There will, of course, be a family of degenerate scalar harmonics 
with K = 2. Some of these will preserve some of the symmetry of CP^ whereas others break it 
completely (this is proved in Appendix B.3[ )p| For a mode of the latter type, the associated metric 
perturbation will break completely the SU{N + 1) subgroup of the R x U{N + 1) isometry group 
of the background metric. It preserves only a R x U{1) subgroup corresponding to time-translation 
invariance and invariance under the rotations generated by d/dip. Hence the corresponding family of 
new black hole solutions will possess just a single rotational symmetry. 

Very recently, Ref. [7j has constructed approximate black ring solutions with just a single rotational 
symmetry. Our results are the first evidence for the existence of new black hole solutions with a single 
rotational symmetry and horizons of spherical topology. 

Note that if one used a k = 2 harmonic that does preserve some of the symmetry of CP^ then 
presumably this would give rise to a different family of new black hole solutions, with more than 
one rotational symmetry. Therefore, assuming that each stationary zero-mode corresponds to a new 
nonlinear stationary black hole solution, there must exist several new black hole solutions that bifurcate 
from the MP family at the same point, and these different solutions have different numbers of rotational 
symmetries. So how many new solutions are there? 

One way of addressing this question is to determine the number of parameters in the most general 
K = 2 harmonic. If we take D = 9 then k = 2 harmonics correspond to the [2,0,2] representation of 
SU{4), which is 84 dimensional. Hence the most general k = 2 harmonic is labelled by 84 parameters. 
Some such harmonics are related by acting with SU (4), i.e. by rotations of the background spacetime. 
However, since SU{A) has dimension 15, this can eliminate only 15 parameters, leaving 84 — 15 = 69 
parameters that cannot be eliminated by rotations of the background. So, up to rotations of the 
background, we have a family of stationary zero-modes with 69 parameters, and presumably a family of 
new black holes with 70 parameters, the extra parameter being the mass (or rjvf )• This is considerably 
more parameters than the 5 that are required to specify the D = 9 MP solution! 



^ We expect that this new family will have unequal angular momenta in general. However, this cannot be seen from 
our results since n = 2 modes do not change the mass or angular momenta at the linearized level (see Appendix E.4I. A 
second order calculation would be required to determine these changes. 

® It is helpful to think about the case of CP^ = , for which k — £, the total angular momentum quantum number. 
Scalar harmonics are labelled by £ and m. Certain modes with £ — 2 (i.e. quadrupole modes) may preserve some 
symmetry (e.g. if m = then they are axisymmetric) but generically they break all of the continuous symmetries of S^. 
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Figure 5: Results in D = 9, k = 2. We represent this mode for fixed values of Fr^ (first graph) 
and ajru (second graph). As opposed to the k = 1 case, the time-dependent mode extends all the 
way to /c = 0. There is not only a new Gregory-Laflamme instability of the black string, but also an 
instability of the black hole, appearing at a = 02, where 02/^4/ = 0.6858 > ai/rM = 0.6300. In the 
first graph, the interrupted vertical line to the right corresponds to extremality. In the second graph, 
the curve shrinks to the origin as a — )• 02. 



4 Thermodynamics and black string instabilities 

4.1 Introduction 

In this Section, we shall explain how many of our results, particularly for k = 0, 1, can be explained on 
the basis of thermodynamic arguments. In particular, we shall explain why the stationary zero- mode 
that appears at a = ai must correspond to a variation of parameters within the MP family. In order 
to do this, we shall develop in some detail the connection between negative modes, thermodynamics, 
and classical stability of black strings. 

We shall consider a rotating black string with metric given by ^ but now with (7^,^ a general MP 
black hole, not necessarily cohomogeneity-1. 

4.2 Rotating black strings 

The black string is expected to suffer from a classical Gregory-Laflamme instability |12j . that is, there 
will exist perturbations which grow exponentially in time as e~^^^ where Im(a;) > 0. This is expected 
to be a long- wavelength instability, i.e. unstable modes have \k\ < kc for some critical wavenumber 
kc- For a static black string, unstable modes have Re(a;) = 0. In the limit k — )• kc, unstable modes 
reduce to a static perturbation of the string, which we shall call the threshold unstable mode. The 
presence of this mode indicates the existence of a new branch of non-uniform black string solutions 
that bifurcates from the branch of uniform solutions [IB]. 

For a rotating black string, we would expect unstable modes to have Re(u;) / in general, and 
so one might not expect to find stationary perturbations with k = kc- Nevertheless, Ref. [21J has 
obtained numerically non-uniform rotating black string solutions that do indeed bifurcate from the 
uniform branch (based on cohomogeneity-1 MP solutions) at a point corresponding to a stationary 
perturbation. Hence stationary perturbations do indeed exist. We believe that the reason for this is 
that unstable modes will have Re(a;) = if they are invariant under the rotational symmetry of the 
black hole predicted by the theorems of Refs. |17|l26j. i.e. the symmetry generated by the Killing field 
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Figure 6: Results in D = 9 for k = 2, in the limit A; = 0. We represent the black hole instability in a 
plot of V tm versus ajru- The interrupted line corresponds to extremality. Numerical error prevents 
us from extending our results all the way to extremality. 



rjjmj where mj are the rotational Killing vector fields and the associated angular velocities of the 
horizon. We do not have a proof of this, but our results, and the results of Refs. |21i B]. indicate that 
it is true. In the limit k — )• fcci this gives a stationary threshold mode that preserves this symmetry. 
This symmetry is a necessary condition for the threshold mode to correspond to a bifurcation into a 
new family of non-uniform rotating black string solutions, since presumably this new family should 
respect the theorems of Refs. |171 126j (although, strictly speaking, these theorems apply only to black 
holes, not black strings). 

In summary, we expect an instability of the black string to appear for wavenumbers \k\ < k^- If 
we restrict attention to modes invariant under the symmetry generated by rijmj then unstable modes 
will have Re(a;) = 0, and the threshold unstable mode, with k = kc, will be stationary and invariant 
under the same symmetry. 

For cohomogeneity-1 black holes, Ojm-j is proportional to d/dtp, so modes invariant under the 
symmetry generated by Jljmj must have m = 0, which is why we set m = above. 

4.3 Thermodynamic negative modes 

Equation ^ shows that a threshold unstable mode of the black string corresponds to a stationary 
negative mode of the Lorentzian Lichnerowicz operator of the black hole. However, in order to make 
contact with thermodynamics, we need to relate it to a negative mode of the Lichnerowicz operator 
for the Euclideanized black hole geometry as was done in Ref. [31] for static black strings. The 
Euclideanized geometry is defined by the analytic continuation t = —it. This produces a complex 
metric since the MP solution is non-static. In coordinates with mj = d/d(j)i, regularity requires that 
we identify (r, 0i) ~ (r, (pi + 2tt) ~ (r + /3, -|- where /? = 1/T is the inverse temperature [9j 

If the Euclidean Lichnerowicz operator in this geometry admits a stationary negative mode, i.e. 
one preserved by the symmetries generated by d/dr, then, after the analytic continuation r = it it 
becomes a stationary negative mode of the Lorentzian black hole, i.e. a threshold unstable mode of 
the black string. 

^ Note that, although the metric is complex, the manifold is real. Real coordinates can be defined by setting 
4>i = <i>i ~ '^^iT, SO the identifications are (t, (j>i) ~ (r, (f>i + 2n) ^ [t + [3, (pi). 
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For static black strings, this correspondence between the threshold unstable mode and a Euclidean 
negative mode of the black hole was used in Ref. |31j to justify the Gubser-Mitra conjecture regarding 
black string stability I16|. In the present context, this conjecture asserts that a black string with 
a non-compact z-direction suffers a classical Gregory-Laflamme type instability if, and only if, it is 
locally thermodynamically unstable. Local thermodynamic stability is the statement that the Hessian 
matrix 

must be positive definite, where S is the entropy of the black hole (entropy per unit length of the 
black string) , which is regarded as a function of the charges = (M, Ji) of the black hole^ The idea 
is that if the Hessian fails to be positive definite then it becomes thermodynamically preferred for the 
mass and/or angular momentum to become non-uniformly distributed along the string. 

The idea of Ref. |31j is to relate local thermodynamic stability of the black string to the existence 
of a stationary negative mode of the Euclidean black hole. This is done by constructing a family of 
off-shell geometries for the Euclidean path integral for which one can show that the Euclidean action 
decreases in a certain direction if the black string is locally thermodynamically unstable. It then 
follows that there must exist a Euclidean stationary negative mode and hence a threshold unstable 
mode of the black string. In the present case, this argument can be phrased as follows. 

Consider a Euclideanized black hole solution B{x) uniquely specified by parameters x". Let T(x), 
Oj(x), M(x), etc. denote the temperature, angular velocities, mass, etc. of this solution. We can con- 
struct an off-shell generalization 13{x,y), specified by parameters y° as follows [3] (see also Ref. p7]). 
Assume the same isometrics as B{x), i.e. Killing fields d/dr and d/d(j)i. Perform an ADM decomposi- 
tion of the metric, with time coordinate r. Take the spatial geometry of B{x, y) to be the same as that 
of B{y). Now choose the lapse function and shift vector so that (i) B{x,y) has the same asymptotics 
as B{x), (ii) B{x,y) is regular everywhere, in particular at the bolt (Euclidean horizon), subject to 
the identifications (r, (pi) ~ (r, + 2-k) ~ (r + /3(x), pi + \ Q-i{x)j3{x)); (iii) B{x^ x) = B{x). Note that 
(ii) implies that B{x, y) is a configuration in the Euclidean path integral defined for temperature T{x) 
and angular velocities J7i(x), for which the saddle point is B{x). Calculating the Euclidean action of 
B{x,y) using the Hamiltonian formalism gives 

I{x, y) = P{x)M{y) - S{y) - p{x)ni{x)Ji{y). (12) 

Condition (iii) implies that the geometry with y = x satisfies the equations of motion and hence the 
first derivative of the action with respect to y" must vanish for y° = x". This is a consequence of the 
fact that the black hole satisfies the first law of thermodynamics, 

dM = TdS + iliJi. (13) 

The second derivative of the action, i.e. the Hessian of the action, now reduces to 

d'l \ ^ ( d'M _ d'S _ d'J. \ 
dy'^dy^)y=x \ dx'^dxP dx'^dx^ ^ 'dx'^dx^)' ^ ' 

where the RHS is evaluated at x. For the MP black hole, the charges M and Ji uniquely parameterize 
the solution (whereas T and Vti do not). Hence we can choose = (M, Ji). We then have 

Therefore, if —Sap fails to be positive definite for some MP black hole then the Euclidean action 
decreases in some direction and hence the black hole must admit a negative mode|^ Given that our 

*This condition can be shown to be equivalent to the standard stability condition in the grand-canonical ensemble: 
the Hessian matrix d^{~G)/dx°'dx^ must be positive definite, where G = M — TS — Q,iJi is the Gibbs free energy 
regarded as a function of x" = (T, Q,i). 

^ Note that we have not constructed the negative mode explicitly by this argument: the linearization of B{x,y) 
around y — x will give a superposition of eigenfunctions of the Lichnerowicz operator. The point is that, since the action 
decreases in some direction, this must involve a negative mode. 
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off-shell geometries are stationary, this negative mode must also be stationary^" In general, there 
must be at least as many negative modes as there are negative eigenvalues of this Hessian. 

We shall refer to a negative mode whose existence is predicted by this thermodynamic argument 
as a thermodynamic negative mode. We shall see below that there are some negative modes whose 
existence cannot be predicted by this thermodynamic argument, so we shall refer to these as non- 
thermodynamic negative modes. Although the former type are relevant for black string stability, it is 
the latter type that are relevant for black hole stability. 



4.4 Rotating black strings are unstable 

For the MP solution, indeed for any vacuum black hole solution specified by (M, Jj), we can now prove 
that —Safs always admits a negative eigenvalue. The proof goes as follows. Consider the Legendre 
transform of the entropy, 

W = S-PM + pQiJi . (16) 



The first law (13) implies that dS/dx'^ = , and also dW/dz"^ = — , where = {M,Ji) and 
= (/3, — The Hessians —Sajs and Wq/j = d'^W/dz'^dz^ are thus inverse matrices, and —Sa/s 
has a negative engenvalue if and only if Wap does. A sufficient condition for Wa/3 not being positive 
definite is 

d'^Wiz'') fdM\ 

The Smarr relation, valid for asymptotically flat vacuum black holes, reads 

^ ~^ M = TS -^iJi ^ M = -{D-2)WT. (18) 



D 



Hence we have 



Woo = -{D -3)MT <0, (19) 



which implies that —Sap always admits a negative eigenvalue. Therefore any vacuum black hole 
solution must admit a stationary negative mode, and a black string based on such a black hole solution 
must always be classically unstable. This explains our result (for cohomogeneity-1 black holes) that 
the black string has an instability in the «; = sector for all values of a. We have argued that a 
stationary negative mode must be present for all values of a, and thus there is no value for a for which 
it reduces to a zero mode. Hence there cannot be an instability of the black hole in this sector. 



4.5 Zero-modes 

We have seen how local thermodynamic instability implies the existence of a negative mode. However, 
we have not yet shown that this negative mode appears precisely when an eigenvalue of —Sap changes 
sign. We shall now address this point by considering zero modes arising from variation of parameters 
with the MP family of solutions. 

Consider, in Euclidean signature, a variation of the parameters x" of the MP solution: x" — >■ 
+ In general, this will induce a variation in T and Oj. Let the perturbed geometry be 

g^iv + h^u- Of course, this is regular provided we make coordinate identifications associated to the 
new values of T and VLi. However, we are demanding that g^y and h^y must be regular separately. 
Regularity of g^y requires that T and ilj must be those associated to the background geometry. But 
then g^j^jj + /i^i, cannot be regular, i.e. h^^ cannot be regular. Hence any variation 5x°^ that leads to 
a change (to first order) in T and must necessarily result in a h^y that is not regular. Therefore 
a variation in parameters in the MP solution gives a regular Euclidean zero-mode if, and only if, it 
preserves T and ilj to first order. 

^° The operator commutes with d/dr so one can work with simultaneous eigenfunctions of these operators. 
Eigenfunctions with different eigenvalues of the latter will be orthogonal. 
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If the variation preserves T and Qi then (using the first law) it preserves [dS / dM)j = 1/T and 
{dS/dJi)M = —^i/T, i.e. it preserves dS/dx°'. Hence it must be a eigenvector of Sai3 with eigenvalue 
zero: 

= 5{daS) = 5x^dpdc,S = S^pSx^ ■ (20) 

Clearly this argument can be run backwards, so Sa/s admits an eigenvector with eigenvalue zero if, 
and only if, there exists a change in the black hole parameters that preserves (to first order) T and 

Consider a 1-parameter subfamily of MP black hole solutions labelled by a parameter A such that 
some eigenvalue of — Sq,^ changes from positive to negative as A increases through Aq, with the signs 
of the other eigenvalues fixed. The above arguments shows that the solution with A = Aq admits a 
stationary zero-mode arising from some variation of the parameters of the MP solution that preserves 
T and fij, and solutions with A > Aq must admit a stationary negative mode. The obvious conclusion 
is that the negative mode is continuously connected to the zero-mode, i.e. as A increases, a stationary 
zero- mode appears at A = Aq and this becomes a stationary negative mode for A > Aq. 

In summary, new negative modes emerge from MP solutions for which an eigenvalue of —Sai3 
changes from positive to negative. Hence a new Gregory-Laflamme instability of the associated black 
string emerges whenever such a change in the sign of an eigenvalue occurs. This is a refinement of the 
Gubser-Mitra conjecture: we are asserting not just that local thermodynamic instability implies the 
existence of a classical instability of the string, but that there is a distinct classical instability present 
for each negative eigenvalue of —Sap- 



4.6 Ultraspinning black holes 

Next we investigate whether —Sap can admit more than one negative eigenvalue. Note that we can 
relate —Sap to the Hessian Hij = {d'^{—S)/dJidJj)M = —Sij defined in Ref. [4J using the following 
identity: 

det{-Sap) = - (^_3)^r dem^) , (21) 

valid for asymptotically flat vacuum black holes. This is proved in Appendix [A} It follows that, for a 
black hole parameterized by (M, Jj), additional negative eigenvalues of —Sap correspond precisely to 
negative eigenvalues of Hij. 

For a MP black hole, for fixed M, the eigenvalues of Hij are all positive for small enough angular 
momenta. However, as some or all of the angular momenta are increased, an eigenvalue of Hij may 
become negative. If we consider the space parameterized by Jj (at fixed M), there is some region 
containing the origin in which Hij is positive definite. We define the boundary of this region to be the 
ultraspinning surface. Following Ref. [3] , we shall say that a given black hole is ultraspinning if it lies 
outside the ultraspinning surface, i.e. if Hij is not positive definite. From the above arguments, we 
know that as one crosses the ultraspinning surface, the black hole will develop a new negative mode. 
The associated black string will develop a new classical instability as this surface is crossed. This is in 
addition to the instability already present at low angular momenta. Furthermore, on the ultraspinning 
surface, the new negative mode must reduce to a stationary zero-mode that corresponds to a variation 
of parameters within the MP family of solutions. 

Let us now examine the form of Hij for MP solutions with different D. First D = 5. In this case, 
there is no ultraspinning region [4J , i.e. Hij is always positive definite. For D = 6,7, we find more 
interesting behaviour shown in Figure [7j We parameterize the black hole by the horizon radius r_|_ 
and the spin parameters ~ Ji/M of Ref. [3D] (rather than M and Ji) in order to produce clearer 
figures. For D = 6, Hij is a 2 x 2 matrix. In this case, the ultraspinning "surface" is a closed curve 
that encloses the origin. Inside this curve, both eigenvalues of Hij are positive but outside it one is 
positive and one is negative. Note that all extreme black holes are ultraspinning. 

For D = 7, Hij is a 3 x 3 matrix. The ultraspinning surface encloses the origin. However, now 
there is another surface on which a second eigenvalue of Hij change from positive to negative. This 
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Figure 7: Parameter space for MP black holes with D = 6 (left) and D = 7 (right). The black hole 
is labelled by the horizon radius r_|_ and the spin parameters which we assume to be positive. For 
D = 6, the blue curve corresponds to extreme black holes. In the red region, both eigenvalues of Hij 
are positive. In the blue region, corresponding to ultraspinning black holes, one eigenvalue is positive 
and the other is negative. For D = 7, the blue surface corresponds to extreme black holes. The 
ultraspinning surface is the red surface near the origin. Inside this surface, Hij is positive definite. 
The orange surface is where another eigenvalue of Hij vanishes. Between the red and orange surfaces, 
two eigenvalues of Hij are positive and one is negative. Between the orange and blue surfaces, one 
eigenvalue of Hij is positive and two are negative. Ultraspinning black holes correspond to points 
between the red and blue surfaces. Note that the "cusp" where the red and orange surfaces meet has 
equal a,, i.e. it corresponds to a cohomogeneity-1 black hole. 



surface lies between the ultraspinning surface and the surface corresponding to extreme black holes, 
and touches the ultraspinning surface at a point with equal aj. For generic Oj, if one gradually scales 
up the Ui then one eigenvalue of Hij becomes negative as the ultraspinning surface is crossed, and 
another eigenvalue becomes negative as the other surface is crossed. As each surface is crossed, the 
black hole should develop a new stationary negative mode, and the associated black string should 
develop a new instability. 

Some D = 7 black holes are special, e.g. the second surface does not intersect the region of 
parameter space corresponding to black holes for which one of the angular momenta (i.e. one of the 
Oj's) vanishes. For example, consider the singly spinning case. All such black holes must possess the 



thermodynamic negative mode predicted by Section 4.4, a second negative mode should appear at 
the ultraspinning surface, and there should be no further thermodynamic negative modes. This is 
consistent with the results of Ref. jl], who found that a stationary negative mode does indeed appear 
at the ultraspinning surface. They also found evidence that there are infinitely many stationary non- 
thermodynamic negative modes that appear at larger angular momenta, the first of which corresponds 
to the threshold of instability of the black hole. Note that we have justified the assertion of Ref. [3] 
that the zero-mode appearing at the ultraspinning surface does indeed correspond to a variation of 
parameters within the MP family. 

4.7 Cohomogeneity-1 black holes 

Another special case are the cohomogeneity-1 MP black holes studied in this paper, i.e. equal angular 
momenta with D = 2A^ + 3. In this case, Hij has a single eigenvalue that is always positive, associated 
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to the eigenvector (1, 1, . . . , 1), and N degenerate eigenvalues (with eigenvectors related by symmetries 
of the background) which change from positive to negative as the ultraspinning surface is crossed, 
with no further changes of sign at larger angular momenta. Note that only the positive eigenvalue 
corresponds to a variation which preserves the equality of the angular momenta. From the above 
arguments, we expect that new thermodynamic negative modes will emerge only at the unique value 
of the angular momentum corresponding to the ultraspinning surface, there will be precisely of 
these, and they will break some of the symmetries of the background. 

The ultraspinning surface corresponds to a = ai (with vm = 1) and we found in Section [s] that 
new stationary (F = 0) negative modes do indeed emerge at this point. Note that our numerical 
calculations used Lorentzian signature, whereas the above discussion concerned Euclidean negative 



modes. However, in Appendix E.3 we show that a regular Lorentzian negative mode (i.e. a solution 
of equations ([S]) and ([9])) that is stationary and has m = does indeed analytically continue to a 
regular stationary Euclidean negative mode. 

These negative modes correspond to k = 1 harmonics. Since these are thermodynamic negative 
modes, we know that the zero-mode at a = ai must be simply a variation of parameters within the 
MP family, as was asserted in Section [3j 

To see that there are precisely N negative modes emerging at a = ai, we use the fact that k, = 1 
harmonics are in 1-1 correspondence with Killing vector fields on CP^ , so there are (A^ + 1)^ — 1 such 
harmonics (see Appendix B.2). However, some of these are related by rotations of CP^, so we need 



to determine how many parameters can be eliminated by rotations. The counting is the same as for 
SU{N + 1) gauge theory with an adjoint Higgs field. Generically this breaks SU{N + 1) to U{1)'^ 
so we are left with N parameters Hence there are independent negative modes that emerge at 
a = ai, in agreement with the prediction from thermodynamics. 

Our numerical results confirm that the stationary negative mode that emerges at a = ai does indeed 
correspond to the onset of a new instability of the black string in the k = 1 sector, in agreement with 
the refinement of the Gubser-Mitra conjecture discussed above. 

For D > 5, comogeneity-1 black holes are ultraspinning for a > ai and hence might exhibit an 
instability of the form anticipated in Ref. [8]. Moreover, as explained at the end of section [2| as we 
increase D, there is more "space" between the ultraspinning surface and the surface of extremality 
for such black holes. Therefore the likelihood of an instability might be expected to increase with D. 
This is in agreement with our numerical results, which show no sign of any instability for D = 7 but 
confirm that an instability is present for D = 9. Note that this instability does indeed occur inside 
the ultraspinning region. 

The onset of instability is associated to the appearance of a new stationary zero-mode (at a = 02). 
This cannot correspond to a variation of parameters within the MP family (we saw above that this 



is possible only at a = ai). Furthermore, we prove in Appendix E.2 that it cannot be a pure gauge 
mode. This zero-mode is continuously connected to a stationary (F = 0) negative mode that exists for 
a > a2. This is an example of a non-thermodynamic negative mode, i.e. one which is not associated 
to an eigenvalue of —Sap being negative. The same behaviour was observed in Ref. [1], i.e. the onset 
of a classical instability of the black hole is associated to the appearance of a new stationary negative 
mode. 

Ref. [3] found that further non-thermodynamic negative modes appear as the spin of the black hole 
is increased still further. For us, extremality imposes an upper bound on the spin of the black hole 
and we do not find any further negative modes beyond the ones associated to the instability in the 
K = 2 sector. However, we believe that, for larger L>, as well as an instability in the k = 2 sector there 
will be further negative modes in sectors with larger k. These new negative modes will be associated 
to new instabilities of the black hole in sectors with larger k. The stationary zero-modes associated 
to the onset of these instabilities will provide evidence for the bifurcation of new families of black hole 
solutions, involving a large number of parameters, and generically with just one rotational symmetry. 

^^Thanks to David Tong for this argument. 
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5 Scalar perturbations and CP^ harmonics 
5.1 Introduction 

The rest of this paper is devoted to explaining the technical details of our work. We shall start by 
explaining the decomposition of metric perturbations into harmonics on CP^ . 

Metric perturbations can be decomposed into scalar, vector and tensor types according to how 
they transform under isometries of CP^ . Pertubations of different type must decouple from each 
other. The decomposition is explained (for a different problem) in Ref. [25]. Tensor perturbations are 
the simplest, these were discussed in Ref. |24j . We are interested in scalar perturbations, for which 
the perturbation can be expanded in scalar harmonics on CP^ . 



We shall assume that our perturbation has been Fourier decomposed as in equation (10). All of 
our numerical results assume m = 0. However, for the sake of completeness, we shall derive equations 
that are valid for non-zero m. In order to do this, we must address a subtlety (already encountered 
in Ref. I24j). that such a perturbation couples with charge m to the 1-form Aa on CP^ defined in 
Section [2j Hence we must consider charged scalar harmonics on CP^ . First we shall describe these 
harmonics and then explain how to construct gravitational perturbations from them. 

5.2 Charged scalar-derived harmonics in CP^ 

We define the gauge-covariant derivative acting on a charge-m tensor field on CP^ as 

Va = Va-imAa, (22) 
where V is the metric covariant derivative on CP^ . 

5.2.1 Scalars 

Charged scalar fields on CP^ can be expanded in terms of charged scalar harmonics defined by 

(p2 + A)Y = . (23) 

These were studied in [19] (see summary in Appendix [B]) , where it is found that 

\ = i{l + 2N)-m^ , i = 2K + \m\ , (24) 

with K = 0,1,2,.... The modulus sign guarantees that the eigenvalue is the same for positive and 
negative charges, the corresponding eigenfunctions being related by complex conjugation. 
Notice that the presence of the 'gauge field' A leads to 

[Va,Vi,]Y = -i2mJabY . (25) 

5.2.2 Scalar-derived 1-forms 



Given a scalar harmonic Y, we can define 



Y„ = --^PaY, (26) 

which transforms as a charged 1-form on CP^ . This can be decomposed into its (1, 0) and (0, 1) parts 
using the complex structure on CP^ . Denote these as and Y^ respectively, where 

J„^Y± = TiY±. (27) 



'^^Note that A = if, and only if, k = m = 0, for which Y is uncharged and constant. In this case there are no 
scalar-derived vectors nor scalar-derived tensors. 
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We shall refer to as scalar-derived 1-form harmonics. We find that they satisfy 

p2Y± = - [A - 2(iV + 1) ^ 4 m] (28) 

and 

OT± = -i=(AT2miV) Y. (29) 

We shall make use of the result that Killing vectors of CP^ are in one-to-one correspondence with 
uncharged (m = 0) scalar harmonics with k = 1 (see e.g. |19] and our Appendix [B|). Given such a 
harmonic Y, the corresponding Killing vector field is — i(Y+ — Y~). 



5.2.3 Scalar-derived tensors 

Following Ref. [25], we decompose a symmetric tensor Yab into its Hermitian (or (1,1)) and anti- 
Hermitian components according to the eigenvalue of the map 

{JYU = Ja'Jb'^^cd . (30) 

If the eigenvalue is +1 the corresponding eigentensor is called Hermitian, and if it is —1 the eigentensor 
is called anti-Hermitian. In the anti-Hermitian case, we can further distinguish between the (2, 0) and 
(0,2) components of Yab, which are defined by J^^Ycb = ^pi'^^'ab with the upper and lower signs for 
the (2, 0) and (0, 2) components respectively. 

The following quantities form a basis for anti-hermitian scalar-derived tensors 
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= K^t) ^ ^ab = ^ia^b) ' (31) 

Y^ denotes the scalar-derived 1-form harmonics of the previous Section, and denotes the projection 
of Da onto its (1,0) and (0, 1) components. Notice that the correspondence between m = 0, At = 1 
scalar harmonics and Killing vector fields on CP^ implies that Y^^^^ vanish for such harmonics. 

Hermitian scalar-derived tensors can be written in terms of a trace, and a traceless part, for which 
the following quantities give a basis 
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9abY, Y+- = V^Y-^ + V^X) - ^ 9ab {V-Y). (32) 
These tensor harmonics satisfy 

= - [A - 4(iV + 3) T 8 m] Y±± , 



^'Y+- = -(A-4iV)Y+-, 

with 

D'Y- = _>--'(A^ + l)T2m(iV + 2) 

C(l 2 CL ' 

-^'^ta = [(A + 2 m iV)Y+ + (A - 2 m N)Ya] 



(33) 



(34) 



5.3 Decomposition of perturbations in scalar-derived harmonics 

Let us consider now the perturbations of the full spacetime metric. We introduce the orthonormal 
basis 

e(0) = fdt, e(^) =gdr, e^^) = h {dil) + A - ^dt) , e^^ = r e^^ , (35) 



^'^This follows from the scalar part of equation (39) of Ref. |25j . 

^'^ In Ref. [2S], hermitian tensors were converted into (1, l)-forms by contracting with J";,. The two quantities written 
here correspond to terms of the form JY and (the primitive part of) dcTY in equation (47) of Ref. [25j . 
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where e*^'-* is the tetrad of the CP^ manifold. The dual basis is then 
1 111 

6(0) = ^{dt + ^ d^) , 6(1) = - (9r , 6(2) = 7- ' ^{i) = " [^(i) " ^(i)) 5^] . (36) 

/ gar 

Take e^"^^ = {e^^^ e^^-*, e^^^} and a coordinate basis dx" on CP^ . The components /iab of the metric 
perturbation transform as scalars under isometrics of CP^ and can therefore be decomposed using 
scalar harmonics on CP^ . Similarly, since we are restricting attention to scalar-type perturbations, 
components of the form h^a and hab can be decomposed using scalar-derived 1-forms and scalar-derived 
tensors on CP^ : 

hAB = Iab Y , 
hAa = r{fXYt + nY-) , 



(37) 

— - ^ (^++ Y++ + H- Y-- + H+- Y+-) + ffa.Y , 

where = {VF^, X^, Z^}, and the functions multiplying the harmonics depend only on (t, r, ^/^) and 
not on the coordinates of CP^ . The real spacetime metric perturbation is given by Re {h^y). Since 
dt and are Killing vectors of the background solution, we will Fourier expand all of these functions 
in t and ip, i.e. we assume a dependence 6"^^*+'™^. It remains to determine the dependence of these 
functions on r. The stability problem will thus be reduced to a system of linear ordinary differential 
equations. We give these differential equations in Appendices [C] and |Dj 

5.4 Boundary conditions 

The metric perturbations must be regular on the future event horizon H'^ . This boundary condition 
can be imposed by considering a basis which is regular on "H^, since the components of the perturbation 
in that basis must be regular. Let us change to the ingoing Eddington-Finkelstein coordinates that 
are regular at 7^+: 

dt ^ dv — y dr , dip ^ dip -j- dr , (38) 

and consider the basis {dv, dr, dip + A — Vtdv, dx"}. Denote the components of the metric perturbation 
with respect to this new basis with a bar (e.g. /gg, W^). Our boundary condition is that these 
components should be smooth functions of (f , r, 0, x°') at the horizon 

For this class of black holes, the horizon is located at the largest real root r = r_|_ of A = g{r)^'^. For 
a non-extreme black hole, near the horizon, A(r) = A'(r_|_)(r — r+) -|- 0[(r — r+)^], with A'(r4.) > 0. 
Using the relation /(r) = r/{g{r) h(r)), we find that, near the horizon, the metric components in the 
original basis are related to the components in the new basis by 



2 A// — \ ' •'01 ~ /oo — 

r1A'(r+j r — r+ r+ 



/oo-2/oi + /ii«A'(r+)/n (r-r+), 

^''^V~mT)7^' /l2-/02^^^A2^A^, (39) 

/22 = 7-7 To /22 ) = TT ^ ) 

r+v/A'(r+) ^r-r+ r+ 



^^In other words, we demand that the tensor field should be smooth at T-L^ . This is stronger than the statement 
that the metric perturbation should be regular at the horizon e.g. it excludes the possibility that /i^^ is singular at T-i^ 
in a certain gauge but can be made regular by a gauge transformation. For example, we show in Appendix |E.3| that, in 
the transverse traceless gauge, a perturbation with oj = m = that satisfies our boundary condition cannot change the 
temperature or angular velocity of the black hole. It follows that a variation in the parameters of the MP solution that 
does change T or Q.h will not give a perturbation /i^i, that is smooth at the horizon in this gauge. 
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The functions H~^~^ , H , , Hi associated to the components of the metric perturbation on CP^ 
are the same in the two bases. 

Since the components in the new basis should be smooth at the horizon, the above expressions 
give us boundary conditions on the behaviour of the components in the old basis. In imposing these 
boundary conditions, it is important to remember that, near the horizon, 

— — — (40) 

where a = /i(r+)/(r+A'(r+)) is positive (for non-extreme black holes). Hence, for example, the radial 
dependence of /oo near the horizon must be 

/oo « y-^) (41) 

where F{r) is smooth at r = r_|_. 

When solving numerically the stability equations, it will be necessary to work with the combinations 
that maximize the information on the boundary conditions. For instance, one should work with /02 and 
/12 — /o2) instead of considering only the leading behaviour of /02 and /12, otherwise the information 
that /i2(r+) - /o2('^+) = is lost. 

As for the behaviour of the perturbations at spatial infinity r — )• 00, we are interested in boundary 
conditions that preserve the asymptotic flatness of the spacetime. For perturbations of the black 
string, the equations of motion ^ then imply that all the functions vanish exponentially for large r. 

6 The eigenvalue problem 

The ansatz for the metric perturbation h^^ is given by Eq. (37). We list in Appendix [Pj the compo- 
nents of the Lichnerowicz eigenvalue equation ^ in the tetrad basis ( 35 ) . These consist of sixteen 
coupled second order ordinary differential equations, each one being second order only in one of the 



perturbation functions, e.g. (72a) is second order in /oo. However, six of these functions can be solved 
for in terms of the ten remaining functions and their first derivatives when we impose the transverse 
traceless (TT) gauge conditions, listed in Appendix [Cj 

Notice that the TT conditions completely fix the gauge in Eq. ([9]) when /c > 0, since the action 
of the Lichnerowicz operator on a gauge mode is trivial, AiV(^^^) = 0. As for A; = 0, there are two 
distinct cases. The first is the limit A; — )■ for which F — )■ 0, as can be seen on the right plots of Figs. [T| 
[2| |3] and |4} The limiting perturbation = is an unphysical pure gauge mode, as happens in the 
original Gregory-Lafiamme case |13j . The second is the much more interesting stationary perturbation 
/cc = marking the onset of a new Gregory-Lafiamme instability when the rotation increases. In the 
left plots of Figs.[3j|4]and[5} this is the threshold mode of the curve F = 0. In the right plots of the same 
Figures, it would correspond to the graph squeezing into the origin for a critical value of the rotation. 
These stationary perturbations are physical since there is no gauge ambiguity. The TT conditions 
require that any gauge vector is a harmonic 1-form, satisfying V^^^ = and VVp^^ = 0. In 



Appendix E.2, we show that no regular harmonic 1-forms exist in this case: they lead to pure gauge 
metric perturbations that diverge either at the boundary r = r_|_ or at infinity r — )■ 00. A proof along 
the same lines, but much more cumbersome, can be given for the modes which represent the actual 
instability of the black holes, i.e. the exponential growth of the perturbations with time (A: = and 
F > 0). 



The tetrad basis (35) is very convenient in the explicit derivation of the TT gauge conditions 
and the Lichnerowicz equations. However, for the actual implementation of the numerical problem, 
it is convenient to choose perturbation functions that make the final equations more amenable to 
numerics, e.g. it is useful to avoid using expressions involving square roots. It is also helpful to define 
combinations of the original perturbation functions which can be solved for algebraically through the 
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TT gauge conditions. Both features are respected if we consider the perturbations in a related basis 
such that: 



foo 
fii 
fo = 
fi = 

h = 
p = 

u = 



f00f-2f02fhn + f22h^n\ hl = foifg- fughn, h2 = f02fh- f22fn, 
fll g"^ , fl2 = fugh, f22 = f22 h'^ , 

r {{w+ + W-) f-{z+ + z-)hn) , fo = -i i r {{w+ - W') f - {z+ - z-)hn) 

-^rg{X+ + X-), h = -ilrgiX+ -X-), 



-\rh{Z+ + Z-), f2 
\{2H+-~H++~H-), 
U2H+- + H++ + H~~) , 



-i\rh{Z+ - Z-), 
Q = i\ {H++ - H- 
V = Hl + ^ {H+ 



(42) 

Now, we solve for six of these functions (foo, fo, f2i Q, U, V) in terms of the ten remaining functions 
and their first derivatives by imposing the TT gauge conditions in Appendix [C| Upon this substitution, 
the corresponding second order Lichnerowicz equations will become third order. The ten equations 
which are second order in foi, fo2, fii, fi2) f22j fo, fi, fi, -f> '^i^ remain second order. They will 
constitute the system of equations to be solved numericallj A non-trivial consistency check on the 
gauge choice procedure is that the ten final second order equations must solve the six third order 
equations (e.g. a third order equation is a derivative of a second order one). We verified explicitly 
that this is the case. 

The final system will be solved using a spectral numerical method. The application of the method 
is simpler for Dirichlet boundary conditions. We consider then the following perturbation functions: 



qi = 

Q3 = 
q5 = 
Q7 = 
qs = 

= 

qio 
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q2 



r J 



1 



r / 



foi J ^ 

fi, c 

p, 

f02 + f22 + 
1 



\ mf2jf )+l 
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fl, 



(43) 



fl2 + 



H 



a Qh 



r 



01 



a 



r 



foi 
fll 



a 



r 



rii 



o + VLHh- — [l 
a 



H]. 



which vanish linearly at the horizon location r = r+. This behaviour can be verified in the expressions 
(39) and (40). The particular combinations chosen for gg; 99 and gio encode the total information 
about the boundary conditions imposed by regularity, as argued in the end of Section 5^ For the 
numerical implementation, it is convenient to use the variable 



y 



1 



r 



(44) 



instead of the radial coordinate r, since y is dimensionless and bounded, < y < 1 . The functions 
represented above vanish at infinity r = oo (y = 1) since the large r behaviour of the solutions of ^ is 

^®The explicit representation of the system would be far too cumbersome. The interested reader can obtain it from 
the TT conditions in Appendix [cj the equations of Appendix [P] and the relations (42 I. 
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exponential, e^^^ , and regularity at infinity imposes the minus sign. The A; = case will be obtained 
as the limit A; — )> 0. 

The system often second order ordinary differential equations is ready to be solved numerically. We 
now describe the numerical method that we employed to solve it. Amongst the most popular choices, 
finite difference methods, finite element methods and spectral methods, we decided to use spectral 
methods. These methods are particularly useful to solve a system of coupled ordinary differential 
equations (ODEs), or coupled partial differential equations (PDEs), to high accuracy on a simple 
finite domain, as long as the data defining the problem are smooth. To our knowledge, the first 
application of these methods to general relativity was given in \28\ H] . 

Typically, in spectral methods one tries to approximate a given function, defined on a finite do- 
main, as a finite sum of algebraic polynomials p{y) = Y^'iLo'^iV^^ i-^- constructs an interpolation 
polynomial. At first sight, it seems that the best idea is to use polynomial interpolation in an equidis- 
tant grid with M + 1 points. However, it turns out that this is catastrophic, because of the Runge 
phenomenon. Not only this approximation does not converge in general, as — )• oo, but also the 
approximation gets worse at a rate that can be as large as 2^ . The correct approach is to perform 
the interpolation in a non-uniform grid. There are many sets of points that are effective, we used the 
so called Chebyshev points given by 

a + b a — b ( jTr\ , , 

%■ = — + ^cos(^^j, iG{0,l,...,AA}, (45) 

for Uj £ [a,b]. Not only does this set of points avoids the Runge phenomenon, by clustering points 
near the boundary, but it has another well-known advantage over uniform grids: it typically leads 
to methods with exponential accuracy. However, because we are effectively approximating a given 
function as a sum of polynomials, we must require our functions to be analytic. 

The procedure to solve differential equations is in the same spirit as standard quantum mechanics. 
Consider a complicated system of n coupled linear ODEs with variable coefficients, 

n n 

J^Hafsqf = -k^^T^pqf\ a E {0, 1, . . . , n}, (46) 

/3=1 13=1 

where each Ha/^ is a second order operator in y, each Ta/s is a scalar function and {k,qjj^} are the 
eigenvalues and eigenfunctions that we want to determine. We now perform our approximation and 



discretise the [a, b] interval in a set of discrete points given by Eq. (45 ). Each q^g^ is then approximated 



(k) 

by a vector, qj^ , whose entries are the values of the eigenfunctions we want to determine at yj. 
Following this procedure, one represents derivatives with respect to y by matrices, Dj\f, that act on 
the vectors q^i^\ mixing adjacent points (see p. 53 of [36] for an explicit construction of such matrices). 
After this approach is complete, each Has and Tq,« are transformed into square matrices H_ and 

X , respectively, of dimension (AA + 1) x (A/" +1). 

All we are left with are the boundary conditions. However we conveniently choose to work with 



the functions defined in Eq. (43), in terms of which the boundary conditions are just of the Dirichlet 



type. One can easily impose this type of boundary conditions by setting the first and last element 

(k) 

of each q^ to zero, which amounts to delete the first and last elements of each eigenvector, and to 



eliminate the first and last column and row of each matrix H and T [36j . In the end we are left 

a/3 af5 

with the following system of linear algebraic equations 

n n 

where I£ and X are obtained from J£ and X by deleting their first and last column and row, 

a/3 0/3 o/3 a/3 



and thus are square (AA — 1) x (A" — 1) matrices. Finally, the system of equations Eq. (47) can be 
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recast in a more convenient form, 



H ... H 

11 In 



K ■■■ K 

nl nn 





- - 











T ... T 

11 In 



nl nn 





- - 















(48) 



which is just a standard generahzed eigenvalue problem in —k'^ of dimension n{M — 1). 

The success of the application of spectral allocation methods to our system of ODEs hinges on 
two key properties that our system of ODEs possess, namely ellipticity of the second order operator, 
meaning that we are dealing with a boundary value problem, and the fact that our system reduces to 
a generalized eigenvalue problem in —k"^. 
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A Determinants of the thermodynamic Hessians 

We wish to show that 

det(-S„^) = - (^_3)^^ det(i7,,) , (49) 

where S^p = d'^S/dx'^dx^ , with = (M, J,), and Hij = d^{-S)/dJjdJj = -Sij. The first step is to 
use the determinant identity 

det ^ = (A - C-^BkBi) detidj) . (50) 

In our case, the left-hand side is det{— Sap), and Cij = Hij. We have then det{—Sap) = adet{Hij), 
with 

. - -s„ - H,, s..s„, - - (asj/ [-m-), [-m-), ■ 

where we have used the first law of thermodynamics dS = PdM — /3ClidJi. Since 

(dmj)y'_f djj \ /dj^\ _( dj, \ ( dm,) \ .52) 



we get 



The identity 



di_\ ^ fdi3_\ _ fdi_\ fdJA 

dMjj \dM)p^ KdjJ^ydMjp^ 



(54) 
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and the Soi = "Maxwell relation", {d{pVLj)/dM)j = (9/3/5Jj) further imply that 



dM J 



-1 



dl3 



Following the steps taken between Eqs. \17\ and (|19|), we finally obtain 

< 0. 



1 



{D - Z)MT 



(55) 



(56) 



Those steps require that the Smarr relation ( |18| ) is valid, i.e. they apply only to asymptotically flat 
vacuum black holes. 



B The geometry of CP^ 

B.l The Fubini-Study construction 

We review here the Fubini-Study construction of the Einstein-Kahler metric and Kahler potential on 
CP^ |19j^ This construction allows us to generate iteratively the CP^ metric and potential from 
the knowledge of the metric and potential of CP^~^. 

Take the C^"*"^ manifold with complex coordinates Z"^ and flat metric 

dslN+2 = dZ"^ dZA , (57) 

where the index A runs as A = (0,a), with 1 < a < A^. Introduce inhomogeneous coordinates 
_ ^ct^^o ii^Q patch where Z^ ^ , such that 

ZO = e^^|ZO|, = Z° C" = i?jv Z^ZA = r'^, f = 1 + = + R% ■ (58) 

Furthermore, introduce a new set of {N — 1) inhomogeneous coordinates (0 < i < — 1) such that 
^Af = gi*^/2 |^7V| ^ ^i^u^yi ^ith -u"n" = l. (59) 

The flat metric on C^^^ can then be written as 

ds2N+2 = dr +r drjgAr+i > where dQ.2N+i = \dT + A(^i^^) + dT.j^ (60) 

is the metric on the unit sphere 5*^^+^, and dT?^ is the unit CP^ metric. Written in this way we see 
that 52^+^ is a Hopf fibration of 5^ over CP^. In ([6O]), A(^n) is the CP^ Kahler potential. Explicitly, 
the CP^ metric and Kahler potential are given by (Pjv ^ and < ^'at < Att) 

dill = g.,dx^d.' = + {d^N + 2A(^_,))^ + -f^ d^%^, , 

= ^rf% ^"^^^ + 2^(^-1)) ' (61) 
in terms of the Fubini-Study metric, dT,'j^j_^, and Kahler potential, ^(jv-i)j on the unit CP^~^, 

d^%_, = fj^\dv'dv'-f]^Wv'dvf, 

An-1) = ^ i /^ii dv' - dv') , fN-i = l + v'v\ (62) 



^^We use the coordinates {Tijv, ^iv} that are related with the coordinates {5, r} of [19] through the coordinate trans- 
formation sin^ ^ = R%/{1 + R%) and f = *jv/2. 
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By definition, the Kahler form on CP^, Jn = ^dA^ , is covariantly conserved, ^aJ^^^ = 0, and 
satisfies JaJbc = -Qac ■ 

The lesson from this analysis is that starting from the CP^ fields we can construct iteratively the 
(jpN geometry as well as the complex coordinates that define the embedding of CP^ in C^"*"^. 
CP^ is isomorphic to the 2-sphere 5^, its metric and Kahler potential being given by 

dS? = J {de'^ + sin^ edcj)'^) and ^(i) = ^ cos # . (63) 

Examples of the embedding in C^+^ may be elucidative. For CP^, parameterized by (0, (/), i?2, ^2), 
the map is given by 

{Z\Z\Z^) = ^l£== fl,fi2Cos^eH(^^+^),i?,sin^e^i(*-<^)) , (64) 
while for CP^, parameterized by {9, (p, R2, ^'2, -R3, ^3), the map is 

In order to reproduce the results in Section 5^ it is useful to recall that, for CP^ , 

Rabcd = 3ac9bd — daddbc + JacJbd — J adJbc + abJ cd ■ (66) 



B.2 Scalar harmonics and Killing vectors 

We review here the systematic way to construct all scalar harmonics and Killing vector fields on CP^ 
[19] . The isometry group of CP^ is SU{N + l). Let Ta^.-Ap^^'"^" be a constant Hermitean SU{N + 1) 
tensor, which is symmetric in the index set {Ai, . . . , Ap} and the index set {Bi, . . . , .6^}, and traceless 
in any contraction between an Ai and a Bi index. This defines the {p, q) representation of SU (N + 1). 
The charged scalar harmonics are then given by 

Y = r^,...A/^-^' • • • Z^^ Zb,--- Zb, , (67) 



and satisfy the Laplacian (23) for A = 2[2pq + N{p + q)]. We have k = max{p, q} and m = p — q. 
Uncharged scalar harmonics have k = p = q and A = + A^). 

The Killing vectors on an Einstein-Kahler space can be constructed from the uncharged scalar 
harmonics with k = 1 , which have eigenvalue A = 4(1 + N). Indeed all the Killing vectors ^(j) of CP^ 
are generated by the relation 

e» = JfN) dtK=m ' (68) 

where i = 1, . . . , A^(A^ + 2) . That is, setting to zero all but one of the constant components of the 
arbitrary Hermitean traceless tensor we get a k = 1, m = scalar harmonic on CP^ through 



(67). Repeating the exercise for all possible combinations, we generate the N{N -i-2) uncharged scalar 



harmonics Y^^^^.^ , and the associated {N + 1)-^ — 1 = N(N + 2) Killing vectors through (68). There 



are N linearly independent Killing vectors which commute with all the others, thus generating the 
Cartan subgroup U{1)^ of SU{N + 1). 



B.3 Symmetries of k = 2 harmonics 

The symmetry group SU{N + 1) is broken, at least partially, by any linear perturbation hat satisfying 
^^hab 7^ , where ^ is one of the N{N + 2) Killing vectors of CP^ . For scalar type perturbations, if 
7^ for some then the symmetry associated to ^ is broken by the perturbation. 
We explained above how to construct the Killing vectors from the uncharged k = \ scalar har- 
monics. Consider now the most general linear combination of Killing vectors K = '^iCi(,(^i), with 
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i = 1, . . . , N{N + 2) . The entire symmetry group SU{N + 1) is broken by hab if the only solution to 
£/^Y = is Cj = for all i, i.e. K = 0. There are (uncharged) k = 2 harmonics on CP^ {D = 9) for 
which this is true. For reference, we present here a particular example in the coordinate system used 
above: a family of m = 0, k = 2 harmonics with three non-zero continuous parameters, /32 and ^3, 



+/32 



/3i 



1 + cos 61 



RiRi{i + cose) 
2(1 + 



62 
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Rl 



2(l + i?i) 



e|i(0-*2)+^3 



'1 + RIY{1 + Rl) 



(69) 



If one of /3i , /32 , /^s vanishes then this is still a k = 2 harmonic but it preserves some symmetry. 



C Traceless-transverse gauge conditions 

The tetrad components of the TT gauge conditions ([s]) are 

- foo + fll + f22 + 2NHL = 0., 



and 
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(71c) 



(71d) 



(71e) 



D Lichnerowicz eigenvalue equations 

In the following, we list the sixteen non-trivial components of the linearized equations ([9]), 



-k^h 



in the tetrad basis (35), given the ansatz (37) for the metric perturbation. Notice that the TT gauge 



conditions have not been explicitly considered yet. The list of components is: 
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- ^ J, 2 JOO + 1- - /II - „ 2 /02 - 7-^ J22 
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r"^ g N VX 

+ — [Z+{X-2mN) - Z~{X + 2mN)] = -k"^ Hl . (72p) 

E Properties of stationary axisymmetric modes 

In the main body of the paper, we have presented our numerical results for general axisymmetric 
time dependent scalar perturbations. Our numerical code has a continuous limit as a; = iF — >• 0. 
Therefore this particular case was already included in our discussion. However, in this appendix we 
want to have a closer look into stationary axisymmetric modes. The reasons are: (i) we developed an 
independent code for this particular case which confirms the results from our general code with time 
dependence; (ii) we proved that the stationary zero-modes (with a; = and k = kc = 0) are not pure 
gauge modes; (iii) we confirmed that our stationary perturbations preserve the angular velocity and 
temperature of the background geometry; (iv) we found which stationary zero-mode perturbations 
can change the mass and angular momenta of the background solution; finally (v) we want to give 
special attention to the stationary modes, and not just to the time dependent instability, since these 
may indicate bifurcation points to new branches of black hole solutions. 



E.l The sub-sector of stationary perturbations 

As described in Section [6j the stability problem of axisymmetric perturbations with time dependence 
consists of a system of 16 Lichnerowicz eigenvalue equations for 16 unknown functions. Choosing the 
TT gauge reduces the problem to a system of 6 TT gauge conditions and 10 Lichnerowicz equations. 
The procedure is consistent because, as explained in that Section, the latter 10 equations automatically 
imply, through the 6 gauge conditions, that the other 6 Lichnerowicz eigenvalue equations are satisfied. 

When we consider the axisymmetric stationary sub-sector of the perturbations, i.e. m = 0, u = 
ir = we find that the initial system of 16 Lichnerowicz equations decouples into a subsystem of 
10 equations for 10 functions and another subsystem with 6 equations involving only the remaining 
6 perturbation functions. Moreover, we can use the stationarity condition, w = iL = 0, to further 
simplify our system of equations. To see how this is accomplished, let us introduce the harmonics 
associated with the time and azimuthal Killing directions as S = e~*'^*e*'"'^. We can then decompose 
the perturbations according to how they transform under the {t, ip} Killing isometrics. For example, 
the scalar-derived vector perturbations satisfy hj^i ~ Ja db§> for b = t,ip and the index A running over 
the radial and CP^ coordinates (the bar in b denotes that the 1-form basis is really {dt, dip+A — ^dt}). 
These perturbations, hAt and h^^, must then vanish when we set a; = and m = 0. In our dictionary, 
this amounts to say that we expect the following: foi = fi2 = fo = f2 = 0. The original 10 time 
dependent Lichnerowicz equations then imply that fi = when we set a; = 0. Similarly, the original 
6 TT gauge conditions imply that Q = 0. We are then led to the following conditions 

foi = fi2 = fo = f2 = fi = Q = , (73) 



or, using the map (pi): foi = /12 = 0, W- = = X+ , = and H" = 



With (73) the original system reduces to a closed system of 3 TT gauge conditions and 7 Lich- 



nerowicz equations. In this axisymmetric stationary case, the boundary conditions at the horizon (39) 



32 



reduce to 
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(74) 



± 



where we used = W'^ir — r+) and one further has W 



-W'^, X- = X+ and Z- 



It 



is important to emphasize that (73) aheady encode the information that the perturbations are in the 



TT gauge, and that the boundary conditions ( 74 ) are compatible with the TT gauge. 



We have done an exphcit search of the stationary modes using only the subsystem of 3 TT gauge 



conditions and 7 Lichnerowicz equations described above subject to (74). We recover independently 



the same results we obtain when we set a; = in our time dependent code. 



E.2 Stationary zero-modes are not pure gauge 

As discussed in Section |6j perturbations with > in the TT gauge have all the gauge freedom fixed. 
However, TT perturbations with k = ^ kc and F = are pure gauge modes [l3j. In this subsection, 
we want to confirm that our stationary axisymmetric zero-modes with k = kc = cannot be pure 
gauge modes. Given that for any residual gauge freedom the gauge parameter would be constrained 
to be a harmonic 1-form, we will prove that there is no regular harmonic 1-form that could generate 
our perturbations. 

Consider the effect of a scalar gauge transformation on the metric perturbations. The most general 
scalar-type gauge parameter can be decomposed as 
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Under a gauge transformation, 

h^u — > hfj_i, + 2 V(^^i,) , 
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(75) 
(76) 



(77) 
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Our stationary axisymmetric perturbations must satisfy (73) which requires that a potentially dan- 
gerous gauge parameter ^ must satisfy 



Co (r) = 6 (r) = , and r{r)=t (r) ■ 



(78) 



We now prove that a parameter ^ obeying these conditions cannot generate a pure gauge metric 
perturbation that is regular. By regularity we mean that the gauge transformation cannot diverge at 
the horizon r = r_|_ nor at the asymptotic boundary r — )• cxd. 

A TT gauge perturbation generated by ^ must satisfy the conditions V^^^ = and D^i, = 0. If we 
introduce the antisymmetric tensor F^j^ = Vj^^j^], for a Ricci flat background, these conditions reduce 

which read simply 



to v^e 



Using 



and V^Pt"" 



^F^"") = . 



„2N+l 



"v/^ and Eq. ( 23 ) , the first of the equations above requires that 



(79) 



(80) 



where the background function g{r) is defined in ([l]) and A is the CP^ eigenvalue (24). The second 
family of equations in (79) further demands that {^~^)' (r) = ^i(r). Introducing the new variable 



6(0 



^(N+3/2)gHr) 



the solution of ( 79 ) must then solve 



x"ii~) = V{r)xi''') ) with V{r) 



Xir) , 



> 0. 



(81) 



(82) 



To study the regularity of the associated gauge transformation, we need the asymptotic behavior 
of x('^) cit the horizon and at infinity. The solution to Eq. (82) can be obtained in these regions by 
considering the dominant contributions of V{r) or by doing a Frobenius analysis. Using (24), we find 
that 



X{r)\ 
X{r)\ 



ao(r-r+) or x(OL~oo, 



(83) 
(84) 



Recall that 6h 



fii/ 



-C^g^y. We have to discard the second possibility in (83) because it would 
generate a dependence /ii|^ ~ (r — r_|_)~^, not compatible with the TT boundary conditions (74). On 
the other hand, we have to discard the solution with the positive sign in (84) because it generates a 
perturbation that grows faster than the unperturbed background metric as r — t- 0. The appropriate 

(85) 



boundary conditions for ( 82 ) are thus 

- ao(r 



xyr) 



We can now complete our proof. Notice that 



< 



x'{rf 



X{r)x{r) 



X{r)x"{r) 



F(r)x(r)2 <0, 



where we used ( 84 ) and ( 82 ) . But these relations can be satisfied only for 

x{f) = =^ ^i(r)=0, and C^{r 



0, 



(86) 



(87) 



which in addition to ( 78 ) proves that our regular zero mode perturbations in the TT gauge cannot be 
pure gauge modes. 
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E.3 Stationary modes preserve the temperature and the angular velocities of the 
solution 



In Section |4j we discussed the connection between the classical instability of the black hole and its 
thermodynamics. This contact was built on the claim that the stationary and axisymmetric modes 
that we study preserve the temperature and the angular velocities of the background solution. Here 
we will prove that this is indeed the case. 

We will first compute the angular velocity and the temperature of the unperturbed background 
solution using standard Euclidean methods. Our strategy then is to check that our TT metric perturba- 
tion h^j,ydx^dx^ is a regular symmetric 2-tensor, when expressed in coordinates where the background 
metric is regular, which confirms that they preserve the angular velocity and temperature. 

We start with the computation of the background angular velocity and temperature. This is done 
performing the standard three steps in the background solution: i) a coordinate transformation to 
coordinates (t, ip) that corotate with the black horizon, ii) a Wick rotation of the time coordinate so 
that we work with the Euclidean solution, and iii) a change to a new radial coordinate that zooms the 
geometry in the near-horizon region. That is, we perform the coordinate transformations 

il) = i)-VlHt, t = -iT, r = r+ + , (88) 

with being the angular velocity of the background solution ([2]). A final coordinate transformation, 

r = 2nTHr, with Th = '-^^^ , (89) 

sets the period of r to be the inverse of the horizon temperature and avoids a conical singularity at the 
horizon. The Euclidean sector of the near horizon region of the background solution ([T]) then reads 

dsl ~ df^ + dp"^ + h{r+f [dTjj + Aadx^f + r\gahdx''dx^ . (90) 

which is a manifestly regular geometry. Indeed the polar coordinate singularity can be removed by 
a coordinate transformation into cartesian coordinates, f = Arctan(y/x) and p = y^x^ + y^. This 
concludes our computation of the angular velocity and temperature of the background black hole. 
To study the regularity of the perturbations, start by introducing the manifestly regular 1-forms, 

E"^ = p^df = xdy — y dx , E'' = p dp = x dx + y dy . (91) 



Consider now our TT metric perturbation h^j^^dx^dx'^ . After using the boundary conditions (74), 
which satisfy the TT gauge conditions, we get 

V dxf" dx" ~ ^^/oiY {p^ df^ + V) + {di^ + Aadx^y + i^^^h^E^ + A^dx'') 

+ ^^hiYiE^^ - ^^'-^^E' {W'a^ + W'-) dx- + 2r+EP (X+ + X") dx'^ 
+ 2 (dV- + Aadx") [Z+ + Z-) dx'' 



+ rl 



dx-dx^ 



(92) 

Notice that the first term on the right-hand side ensures that there is no conical singularity if f has 
the same periodicity as the background metric. Furthermore, the remaining dependence on f and p 
is given by the manifestly regular 1-forms E'^ and E^. Hence, h^ydx'^dx" is a regular two tensor in 
the cartesian coordinates (x, y, ip, x") where the background metric is regular, which confirms that the 
perturbations indeed preserve the angular velocity and temperature. 
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E.4 Modes that can change the mass and angular momenta 

The stationary zero-modes can potentially change the mass and angular momenta of the geometry. 
In this subsection, we find which of our perturbations with a; = 0, m = and k = can change these 
conserved charges. 

The change on the conserved charges associated to a Killing generator ^ introduced by a pertur- 
bation h^i, can be defined via a surface boundary integral as (for TT perturbations) [H [2] 

Qdf^^ 9] = I ea/3M. [i.^'h^'' - h^'^V^e + U^'i^'^C. + ^.e)] dx^ A dx^ . (93) 

The conserved charges of interest are the energy, for ^ = —d/dt, and the angular momenta associated 
with the y{D — l)/2j U{1) Killing vectors = d/d^i. The corresponding changes are denoted by £ 
and J7i, respectively. The surface integral is over a constant time hypersurface at asymptotic infinity, 

To compute the charges of our perturbations, we need the asymptotic behaviour of our solutions. 
This can be obtained from a Frobenius analysis of the Lichnerowicz equations at r — )• oo. We are 
interested in boundary conditions that preserve the asymptotic fiatness of the spacetime, i.e. that 
decay (strictly) faster than the background geometry. We find that the behaviour of the regular 
perturbations is such that they decay at infinity according to 

fn « O {r-'-''^) , h,^0 {r-'^) , « O {r~''^) , fi « O {r-'-"-) , 

fo « O (r~2«) ^ « (^-2k^ ^ p ^ (r-2-2«) . (94) 

We then find the generic behaviour for the changes in the conserved charges, 

J^^ = A, r2-2« + Bi r-2« + Q {r-^-^") , (95) 

where Eq, Ai, Bi are functions of k and r_|_ and, in particular, Ai = Bi = for k = 0. 

We thus conclude that the k = modes are the only ones that decay sufficiently slowly to change 
the mass of the geometry. Moreover these modes cannot change the angular momenta. The only 
modes that change the angular momenta are those with k = 1. All other modes with k > 2 do not 
change the mass nor the angular momenta. We have explicitly done these computations for the D = 7 
case. 
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